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ABSTRACT 

The decay B~ — ► K~K~ir + has been sugested as a test for minimal su- 
q . persymmetric standard model and for supersymmetric models with R-parity 

**^ ! violating couplings, in view of its extreme smallness in the standard model. 

Q_i! We calculate two long distance contributions to this decay, that associated 

o ' with DD and Dir intermediate states and that induced by virtual D, tt 

mesons. The branching ratio due to these contributions is 6 x 10 -12 , which 
is somewhat smaller than the standard model short distance result, leaving 



X 



this decay free for the search of new physics. 



The standard model (SM) of strong and electroweak interactions is 
presently in very good shape. The experimental data agree with SM and 
this continuing success is somewhat paradoxically, a principal factor in the 
intensive search for physics beyond the standard model. This search is con- 
ducted nowadays in various sectors of particles phenomena. Among these, 
rare b decays is considered to provide good opportunities for discovering new 
physics beyond SM 0. 

Among the rare decays studied so far b — > S7 plays prominent role. 
The measured rate in two different experiments is Br(b — > 57) = [3.15 ± 
0.35(stat)± 0.32(syst) ± 0.2Q(mod)] x 1(T 4 § and Br{b -> 57) = [3.11 ± 
0.80(stat)± 0.72(syst)] x 10~ 4 @, to be compared with the latest theoretical 
calculations within the SM giving Br(B -> A s 7) = (3.32 ± 0.30) x 10" 4 [§. 
The agreement with SM is impresive and it is doubtful that a deviation could 
be detected, even when the above figures are improved. This, in view of the 
fact that long-distance (LD) contributions are also present; although more 
difficult to calculate with good accuracy, the existing estimates concur that 
these are approximately (5 — 10)% of the short-distance (SD) amplitude [|]. 

An alternative approach to the identification of virtual effects from new 
particles in b decays like b — > 57 and b — > sl + l~ is the consideration of 
rare decays which are of negligible strength in SM. In such cases the mere 
appearence of the decays at a rate much larger than it is possible in SM 
would be a clear sign of new physics. 

Recently, the decays b — > ssd, b — > dds were proposed [f| as ideal proto- 
types of the latter method. As shown in Ref. 0, the b — > ssd is mediated in 
the SM by box-diagram and its calculation results in a branching ratio nearly 
of 10 -11 , the exact value depending on the relative unknown phase between 
t, c contributions in the box. The b —>■ dds branching ratio is even smaller by 
a factor of about 10 2 , due to the relative |Vtd/Ks| factor in the amplitudes. 
The authors of ref. have also calculated the b —>■ ssd transition in various 
"beyond the SM" models. It appears that for certain plausible values of the 
parameters, this decay may proceed with a branching ratio of 10~ 8 — 10~ 7 
in the minimal supersymmetric standard model and in two Higgs doublet 
models. 

Moreover, when one considers supersymmetric models with R-parity vi- 
olating couplings, it turns out that the existing bounds on the involved cou- 
plings of the superpotential do not provide at present any constraint on the 
b — ► ssd mode 0. It has been pointed out in Ref. || that the hadronic chan- 
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nels most suitable for the search of the b — *> ssd transition are the AS" = 2 
decays B~ — > K~K~n + or B° — > K~ K~7t + tt + . The appropriate exclusive 
channel for b — > dds transition would be B~ — > K+n^n^ . At present there 
is no published experimental limits on these modes. 

In the analysis of Ref. || only the short-distance contributions were con- 
sidered in detail. However, it is well known that long-distance contributions 
which are associated with low-lying intermediate hadronic states are also 
present in particle transitions. As we mentioned above, such contributions 
to B — > Xg'-f are rather small. However, each specific decay mode requires 
the estimation of its LD contribution; this is imperative, since only when 
a trustworthy estimate of such contributions is available one may proceed 
to compare the specific transition to the theoretical SM treatment or use it 
for revealing new physics. This necessity is best exemplified by known oc- 
curences in K - physics M: in some decays like K + — > 7r + 7r°7, K + — > 7r + Z + Z~ 
the SD contribution is obscured by LD contributions while for K° L — > 7r°z/z/, 
K + — > 7i + uu, the LD contributions are considerably smaller than the stan- 
dard model short- distance amplitude ||. The calculation of long-distance 
contributions to a specific process is not based on a well-defined theoretical 
procedure. Clearly, the intermediate states are the main contributions to 
this part of the amplitude. However, the technique of their inclusion, as well 
as the choice of relevant states will influence the final result. We shall rely 
on the accumulated experience from the treatment of long-distance contri- 
butions to various processes, like K — K transition j7|, D — D transition 
[ PHI , K + — > -k + vu [Oj], B — y X s 7 and B s — > 77 [Q decays, in order to 



formulate our approach to the B~ — > K~K~tt + process at hand. 

We include two contributions in the calculation of the long distance am- 
plitude B~ — > K~K~ir + : (I) the box diagram, shown in Fig. 1, which is 
essentially the LD analog of the SD calculation in the standard model 
of the b —> ssd transition. (II) the contribution of virtual " D°" and "7r°" 
mesons, via the chain B~ — > K~" D°" ("71" ") — > K~K~n + . This contribution 
arises as a sequence of two AS = 1 transitions and may lead to final K~K~n + 
state as well. It is therefore necessary to have an estimate of its relevance 
vis - a - vis the "direct" AS = 2 transition. Let us consider firstly the am- 
plitude arising from (I). The two diagrams (a) and (b) express the Glashow 
- Iliopoulos - Maiani symmetry, so that the decay amplitude vanishes in the 
limit m c = m u (m^ = m^). Since each diagram contains two W's, it is re- 
lated to several semileptonic processes with one virtual meson. Thus diagram 
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(a) relates to D° — > K~e + v e , B~ — > D°e~v e and /}° — > 7r + e~i/ e involving the 
product V^V^V^V^* and diagram (b) relates -D — > K~e + v e , B~ — > D°e~v e , 
K~ — > 7r°e~z/ e and 7r + — >■ 7T°e + i/ e involving the product V c bV* s V u dV* s . 
The transition probabilty is given by 

< ^-;r-7r+|s|£r > 6o:r = (-^=)* Vcb v; s v cd v: s J d A qx d%^Q 2 ^Q x 

{5\Qi + q 2 - k 2 )5\ PB -q x - Q 1 )5\Q 2 - Pir - g 2 )5 4 ( gi - Q 2 - A*) 

~ Z( ^~ g2 ^f /M ^ ) [< A-|(gcr|U° > -^-^ < n+D°\(cdr\0 > 
Q2 M w q 2 m D 

- < K-\(suY\7r° > 9 - - < 7i + 7r°\(ud)P\0 >} 



+6\Qi -to- k 2 )5\p B - qi - Q l )5\Q 2 - Pn + q 2 )5\ qi -Q 2 - h) 

< D°\{cby\B- > fi-T < K-\{scY\D» > -^-QiuQiJM^) 
q l — m D — ivi w 



Ql - M w 4-m 2 D 

- < ^-7T°|(^) Q |0 > 2 - . < 7t + \(udf\n° >] + {h <-> k 2 )}, (1) 

where ((jjqi) a stands for q~jj a (l— 75)91, while the rest of the notation is defined 
in Figure 1. The first part comes out from the diagrams on Figure 1, while the 
second results from the crossed diagrams. The calculation of (JJ) depends on 
the matrix elements < D°\(cb) u \B- >, < K'\(sc)^\D° >, < ir + \(cdy\D° >, 
< K~\{suY\ti > and < n°\(ud)P\iT~ >. Since only pseudoscalar states 
appear, we have to deal with transitions between such states induced by the 
vector current only, 

< P'{p')\qn^Pkp) >= UtfW+p" 1 ) + f-{q 2 W -p'n, (2) 

which may be rewritten as ]n| 

<P'(p')\q j l»q i \P(p)> = F l {q 2 ){ P ^ + ^- m2p ~ 2 m2p ' ^-^)) 

+ F (q 2 ) mp (p» - (3) 
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where F% and Fq contain the contribution of vector and scalar states respec- 
tively and q 2 = (p — p') 2 . Also, -Fi(O) = F (0) For these form factors, 
one usually assumes pole dominance JTT], |T3|, [LJ], [T5j 

= F tl f) = ^| (4) 

and in order to simplify, we shall take my = m$, from which results f-(q 2 ) = 
0. We shall assume that one can safely take f+(q 2 ) — 1 pi and the limit 



Qi > Q2 ^ Myy, which then leads to a more tractable expression for the real 
part of the amplitude 

A b r(B-( PB ) - k-^k-^)^^)) = — Vcb v: s v cd v: s j d% 

{ 2 1 2 7 ; ^—ti rl(~ m B + 2k 2 ■ PB){m 2 K + 2A4 ■ jv) 

gf - mj, (g x -k x - p n ) 2 - m 2 D 

+q 2 (m 2 K + 2k x -p^ + ml- 2k 2 ■ p B ) + 2k 2 ■ qi{m 2 K + 2k x ■ p n ) 

+2p w ■ q x 2k 2 ■ q 1 - 2k 2 ■ q x q\ + 2p„ ■ q x q\ + 2p w ■ qi{-m 2 B + 2p v ■ p B ) - qf] 

gf - m£, (gi - A* - p n ) 2 - ml 

+Qi( m K + 2fci ■ Ptt + ^| - 2A; 2 • Pb) + 2/c 2 • + 2k x ■ p w ) 

+2p n ■ q x 2k 2 ■ gi - 2k 2 ■ q x q\ + 2p w ■ q x q\ + 2p n ■ qi(-m 2 B + 2p n ■ p B ) - qf] 

+ (h^k 2 )}, (5) 

where G = v / 2fi ,2 /(8M^ / ). The separate contributions of DD and D% inter- 
madiate states diverge as fourth power. However, the GIM cancellation acts 
in such a way as to decrease the degree of divergence and finally the integral 
in (|5[) will give a quadratic divergence. Similar situations were encountered 
in previous LD calculations [0, [IT| . We note that the explicit inclusion of the 
pole-type form factors for f+(q 2 ) would reduce the degree of divergence, in 
such a case, however, the evaluation of the integrals becomes very cumber- 
some and this effort is not justifed since as it will turn out the contribution 
of the real part is essentially negligible in comparison to that provided by 
the imaginary part. The integrals in ([5]) are calculated by using Feynman 
parametrization. The final result for the decay rate is 

(m B ~m K ) 2 r(si) 



T(B-^K-K-^ = ^-^g^ ds2 j M M (6) 
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where 

(si)i, 2 = m 2 K + ml - ^-[(s 2 -m 2 B + m 2 K )(s 2 + m 2 n - m\) 

(7) 

and A (a, 6, c) = a 2 + b 2 + c 2 — 2a6 — 2ac — 2a6. The A b ° x denotes the lead- 
ing term of the amplitude, which results after using the primitive cut-off 
regular ization: 

Al°%B~ - K-K-1T+) ~ G 2 V cb V; s V cd V; s -±- 2 A 2 (m 2 D -ml). (8) 

There is obviously the uncertainty in the value to be taken for A. The 
momentum in the box cannot exceed ra^, and by taking A ~ 10 GeV we 
obtain 

BR(B~ -> K-K-n + )p } x) ~ 8 x 10~ 15 (9) 
for the real part of this contribution, using T(B~ — > a//) = 4 x 10 -13 GeV 
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Turning now to the imaginary part of the B K K tt + amplitude 
provided by the DD and intermediate states, it is given by 

A h r{B- -> K-K-n+) = -^- Vcb V: s V cd V; s 

62tc z 

J Ai%i - rr&)P((?i - *i - k 2 ) 2 - mf,) - 5(( qi -h- PlT ) 2 - ml)} 

{-qf + -2k 2 ■ qiqf + 2p n ■ q x q\ + q\{m 2 K + 2k\ ■ p n + m B - 2k 2 ■ p B ) 

+2p n ■ q x 2k 2 ■ q 1 + 2p w ■ qi(-m 2 B + 2k 2 ■ p B ) 
+2k 2 ■ qi{m 2 K + 2p„ ■ k{) + (-m| + 2k 2 -p B )(m 2 K + 2k x • p v ) 

+ (h^k 2 )}. (10) 

Introducing now s± = (p B — ki) 2 = (k 2 + p n ) 2 and s 2 = (p B — k 2 ) 2 = 
(ki + P-k) 2 one arrives at 

A b ° x {B- -> i^-7r + ) = --^y^;^ 

x{F(m^, Sl , s 2 ) - F(m 2 , s X) s 2 ) + (s 1 «-> s 2 )}, (11) 
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with 



X 1 l 2 {m 2 D ,m 2 P ,s 1 ) 

F(mp,s 1 ,s 2 ) 



m z D — mp + si 

{-m 4 D + m 2 D (2 Sl - ^m 2 K - ^m 2 - ^s 2 ) + (s 1 - m 7r )(-s 1 + m 2 K ) 

+ (si^s 2 )}. (12) 

Using the expression (|(J) for the decay width we find 

BR(B- -> K-K-7r + )^ ox) - 6x 1(T 12 . (13) 

We proceed now to estimate the second possibility for a LD part which 
may lead to a final K~K~7t + state. This possibility is expressed as two 
consecutive two-body nonleptonic transitions (see Fig. 2) in which the con- 
necting single particles D° and 7r° is virtual. An estimate of this contribution 
requires the knowledge of the < B-\H w \K~ n D°" >, < "_D°" \H w \K~ti + > 
amplitudes for virtual D°, which is lacking. For a virtual ir° existing es- 
timates for < "7T°" \?i w \K~7T + > indicate [TlJ that it is smaller than the 



physical amplitude in a certain region. We rely in our estimation on the 
"physical" amplitudes B~ — > K~D, D° — > K~tt + keeping in mind 



that this induces an amount of uncertainty. However, the final numerical 
results will show that this is of no consequence in the present problem. In 
the diagram (2b) the D° may also be on the mass shell. Therefore, we must 
exclude in our calculations the region around physical D°, which represents 
two AS* = 1 physical decays, B~ — > D°K~ followed by D° — > K~tt + , since 
we are pursuing the B~ — > K~K~tt + outside the resonance region. We shall 
return to this point below. 

The calculation of the virtual D° mediated part of the amplitude requires 
the use of the effective nonleptonic Lagrangian. The part relevant for the 
present calculation is 

C LD = -^{V cb V:^\cbr(su), + a ( 2 b \cunsb),] 

+ V cs V: d [a[ C \csY(du), + aficuYids)^] + h.c.}, (14) 

V qiq2 are CKM matrix elements and a{ \ a 2 , af^ and affl are effective Wilson 
coefficients (see Bauer et al., Ref. |I5|1 ) at the charm and beauty scales. We 
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use factorization approximation [T5j for the two parts of the B — > K K n~ 



amplitude and the expression obtained in [jlj] for the < P%P2\D > transition, 

M <PlP2lD> = ^C PlP2 tfp 2 F D -' p ^m 2 P2 )(m 2 D -m 2 Pl ). (15) 

In (|15D Cp 1 p 2 contains CKM matrix elements and a Wilson coefficient. By 
using the explicit form of ( |T3| ) we have neglected the small contribution from 
the annihilation part of the amplitude, which is proportional to a^, t4 [0 • 
The part of the decay amplitude due to the " D°" pole is then given by 

xF (m K )F (m J q2 _ ^ + —^r • (16) 

The decay width due to this contribution is given by 

r(B- -> K'R-^) = — — | T \C\ 2 \F, BD (m 2 K )F n DK (m 2 )\ 2 

V ; 2(27r) 3 32m|' 110 V KJ V Wl 

(m B -m K ) 2 _ f m | - s) 2 (s - m?^ 2 



x / ds- 



(m^+m K )2 (s - m 2 D ) 2 + (m D r D ) 2 

iA 1 / 2 ( S ,m|,m 2 Y )A 1 / 2 ( S ,m 2 ,m 2 Y ), (17) 
s 

with C = (G 2 /2)V c iy* s V cs V* d af^a^ JkI'tt, for the resonance in a s channel 
and the same for the resonance in a crossed channel. 

Using for af \ af \ Fq K and F BD the values of Bauer, Stech and Wirbel 
[ p~5|j we calculate the virtual contribution by deleting a width of 2 A 

around D mass in the s variable. The size of A is related to the experimental 
accuracy of the D - determination in the final K~tt + state. In the various 
experiments it ranges between 1 and 10 MeV. One should keep in mind that 
average accuracy of D° - mass determination is 0.5 MeV |16|| . Thus, in order 
to delete the physical D 0, s one must take at least A = 1 MeV. However, 
we shall check the A dependence for a range of values to make sure that our 
conclusions are not affected. 

For A = 20, 5, 1, 0.1 MeV we find the nonresonant D° contribution to 



be 

'(DO) 



BR{B- -> K-K-it + ) p , X = (0.31; 1.2; 6.2; 61) x 10~ 15 . (18) 



7 



In all cases the result is much smaller than (13), though one should remember 
that A = (1 — 5) MeV is the realistic option. A similar calculation for tt° 
intermediate contribution, i.e. B~ — > K~ n ix Qn "7r°" — > K~tt + yields a value 
smaller by four orders of magnitude, especially as a result of CKM angles. 
The pole contribution is therefore considerablly smaller than the LD box 
contribution calculated with DD and Dn intermediate states; thus the total 
branching ratio from all diagrams we included is 

BR(B~ -> K-K-n + ) LD = 6 x 1(T 12 . (19) 

As a check, we used our A p ^t\ amplitude to calculate B~ — > K~K~ir + as 
given by decay via a physical D° and we find a branching ratio of 7 x 10~ 6 . 
This agrees very well with the experimental expectation of (9.9 ±2.8) x 10~ 6 , 
obtained by using BR(B~ -> D°K~) = (2.57 ± 0.65 ± 0.32) x 10~ 4 |TJ and 
BR(D° -> i^-7r+) = 3.85 x 10~ 2 [|TJ. 

A few remarks concerning our approximations. As we mentioned, we have 
neglected the form factor dependence in the calculation of A b T ox . Their inclu- 
sion would have decreased the degree of divergence. However, in view of the 
smallness of the real part of the amplitude, this neglect is of no consequence. 
A possibly more serious uncertainty is caused by the fact that we used only 
DD and Dtt intermediate states in the box calculations. Additional interme- 
diate states, all within the physical region could be DD*, D*D*, Dp, D*p, 
Dr](j]'). We did not consider these states for two main reasons: first, there is 
no knowledge of the matrix elements and the required form factors involved. 
Moreover, the inclusion of strongly decaying resonances D*, p as intermedi- 
ate states is questionable and some of their effects are taken into account by 
the form factors considered (|j) . We decided therefore to ignore these contri- 
butions, though we are aware of the possibilty that additional intermediate 
states in Fig. 1 might increase our result by a factor of, say, 2-3. Finally, it is 
interesting to note that our result, which indicates that the LD contribution 
in the B~ — > K~K~n + decay is smaller or at most comparable to the SM 
short-distance contribution, fits into the general picture of B decays. This is 
in contrast to the situation in the strange and charm sectors: in K transi- 
tions the two contributions are comparable in some cases, SD dominates in a 
few decays and LD in many others; in D - transitions LD contributions are 
generally larger than the SD ones, except for the unusual case of B c — > B*^ 
decay @. 
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To summarize, we have shown that the long - distance contributions to 
B~ — > K~K~n + are smaller in the SM than the short - distance box dia- 
gram, and have the branching ratio in a 10~ 12 — 10~ n range. This is a most 
welcome feature since it strengthens the suitability of the B~ — > K~K~n + 
decay as an ideal testing ground for physics beyond the standard model, as 
originally suggested in ref. |J. We expect that this avenue will be explored 
experimentally in the near future and we note that the first analysis of this 
mode has just been completed by the OPAL Collaboration |19| and an upper 
limit of 1.29 x 10 -4 at 90% confidence level has been set for the branching 
ratio of this decay. 
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Figure Captions 

Fig. 1. Long distance box-diagram contributions to B~ — > K~K~n 

Fig. 2. Pole contributions to the long distance amplitude of B 
K~n + , (a) quark picture, (b) hadronic picture. 
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Figure 2: Pole contributions to the long-distance amplitude of B — >■ K K n + . (a) quark picture, 
(b) hadronic picture. 
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